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We study trapping of a cold atom by a single vortex line in an extreme type II superconducting 
chip, allowing for pinning and friction. We evaluate the atom's spin flip rate and its dephasing due 
to the vortex fluctuations in equilibrium and find that they decay rapidly when the distance to the 
vortex exceeds the magnetic penetration length. We find that there are special spin orientations, 
depending on the spin location relative to the vortex, at which spin dephasing is considerably reduced 
while perpendicular directions have a reduced spin flip rate. We also show that the vortex must be 
perpendicular to the surface for a general shape vortex. 
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I. INTRODUCTION 

A significant goal of atom chip experiments is to trap 
cold atoms near a surface at the submicron scale. The 
magnetic fluctuations near a metallic surface lead to sig- 
nificant spin flip (sf ) transition towards untrapped mag- 
netic sublevels and hence loss of atoms from the trapM 
This has led to theoretical study of superconducting atom 
chips predicting a significant reduction of noise^"^ of 6- 
12 orders of magnitude. The reduction is more signif- 
icant when the atom's distance z from the surface is 
in the range A < z < 5 where A is the London pen- 
etration length and S is the skin depth of the normal 
phase; e.g. for Nb chip and Rb atom with sf frequency of 
v — 560KHz we have A = 35nm and i5 = 150/im. Exper- 
imental data^' have reached w 30/im showing, however, 
an enhancement of the lifetime by a factor ss 10. 

In current atom-chip experiments, DC magnetic fields 
of the order of 10— 100G are applied orthogonally to slabs 
of type-II superconductors, resulting in vortices within 
the superconducting material. Dynamics of high density 
vortices were considered as a source of noise^, and the 
relation to flux flow was studiecP. Furthermore, magn etic 
fluctuations lead to dephasing of coherent spin states^!. 

Further interest in vortices is their control of the mag- 
netic field close to the surface hence producing a mag- 
netic trap. Stable traps due to vortices on a thin su- 
perconducting disc were demonstrated experimentally^. 
Other shapes of chip can lead to programmable mag- 
netic trap geometries^. Furthermore, isolated vortices 
can be generated in a remanent state, leading to stable 
traps^l. Near-field noise is expected to be reduced due to 
the proximity of the superconductor, and technical noise 
is minimized as no transport current is needed to create 
the trap. 

In the present work we consider fluctuations of a sin- 
gle vortex and the resulting dephasing and sf rate of an 
atom above the surface. We consider the limit of type II 
superconductors, i.e. the ratio k = A/£ is large, where £ 
is the coherence length. We also assume that dissipation 
is dominated by the vortex frictional motion, i.e. the 
surrounding superconductor is non-dissipative; the latter 



has been studied separateljEH^. We consider the vortex 
displacement u(z), where z < is the distance from the 
surface at z — 0, and show that it must end perpendic- 
ular to the surface for large K, i.e. du/dz\ = 0. While 
this boundary condition has been p reviou sly use d 13 ^ 
it was not explicit in related works^Mil and in fact 
configurations which deviate from this condition were 
considered^. In addition to elastic and external forces, 
the vortex responds to pinning and friction forces. We 
consider strong pinning due to a columnar defect, or 
weak pinning where friction dominates. We find that 
the magnetic fluctuations decay rapidly at distances be- 
yond A, yet even at short distance they have significant 
minima corresponding to a location dependent eigenvec- 
tor. Choosing th e tr ap direction along this eigenvector 
reduces dephasing^l considerably, while choosing it in 
the perpendicular direction reduces the sf rate. 

The strategy is to consider a magnetic dipole m at 
position r = 0, z > (r and u(z) arc 2-dimcnsional vec- 
tors) that emits a magnetic field h l (r, z) with frequency 
u) that is incident on the surface z = 0. The boundary 
conditions determine a reflected magnetic field h r (r, z), 
hence a response function hl(0,Zo) = aijirij. Consider 
a spin polarized in a direction n which is the static mag- 
netic field at the trap center. The fluctuation dissipation 
theorem determines then the magnetic fluctuations, and 
the sf rate of an a tom with moment /i is given by the 
Golden rule asP^ 

2fi 2 J2i Im[ a u] m 
1 s f ff e hu>/k B T _ i ■ W 

Here u) is the transition frequency between the spin levels, 
T is the temperature and i sums on the two perpendic- 
ular directions to n (for spin > 1 one needs 2 or more 
transitions of the form (jTl). In contrast, dephasing is 
caused by the fluctuations in the energy difference of two 
trapped magnetic sublevels^, hence it is determined by 
Im[au] where i in now in the h direction. 

The evaluation of the reflected wave proceeds in the 
following steps: (i) evaluate the magnetic field at z < 
and then h r (r, z) in term of a general vortex shape u(z), 
(ii) find an equation of motion for u(z), including elastic, 
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external source, friction and pinning forces, and solve as 
a response to m, (iii) combine (i) and (ii) to find the 
response a^- and hence the dephasing and sf rates. 



II. MAGNETIC FIELDS 

Consider a superconductor occupying half space at 
z < with a single vortex line whose equilibrium po- 
sition is at r = Uo, i.e. it is perpendicular to the surface. 
Allowing for fluctuations, the vortex position becomes 
u(z) so that the local vortex orientation is z + du/dz and 
the superconducting phase <p(r, z) satisfies 



V x Vip = 2n5 2 (r - u)(z + —) 



(2) 



This equation satisfies § • dl — 2ir with the line inte- 
gral in the r plane around r = u(z) at any given z < 0. 
It is assumed that u(z) is single valued, i.e. the vortex 
does not bend by more than tt/2. Eq. ^ can be solved 



(z+fg)x(r-u) 
(r-u) 2 



(3) 



The vector potential A(r, z) is a solution of London's 
equation^ 



A + A 2 V x V x A 



2tt 



(4) 



where (j>o is the flux quantum. The equation for the mag- 
netic field h(r, z) — V x A(r, z), or its Fourier transform 
h(r, z) = J k h fc (z)e lk ' r where f k = J d 2 k/(2n) 2 , is 



k 2 -^}h k (z) 



V^O -jk-u(z) 

A 2 



.„ du, 

(z+a- z ) (5) 



The solution is found by the Green's function Gk(z,z') 
for h fe (z) = h fc (z) - h fe (0)e'*, with p = 



k 2 , which 



has the desirable boundary condition hfc(0) = 0, 



2p 



z '\ — b p( z + z ')] 



(6) 



The overall solution then is h(r, z) — h p (r, z) + h (r, z) 
where h (r, z) is a solution of the homogenous part of Q, 
to be determined by matching with the external fields. 

We proceed now to study boundary conditions. For 
low frequencies uj the dominant term in the dynamics is 
friction, linear in uj. The Maxwell equation in the vacuum 
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z > is then V x V x A = 0, neglecting the t-A term. 
The current vanishes at the surface, i.e. (V x h) z (r, 0) = 
0, and as usual, h(r, z) is continuous across z — 0. Note 
that V • A = ^V 2 <p9(-z) leads to a jump in d z A z , 
however this boundary condition is not needed for h(r, z) 
(Eq. 10 below) within our quasistatic limit. 

The vector potential at z > has incoming and re- 
flected components 



A(r,z>0)= f A^z)e lk ' r +^ + f 

Jk Jk 



Al(z)j 



k-r— kz 



(8) 



where k = |k|, the gauge is (ik+kz)- A 1 = (ik— kz)-A r = 
and Aj = 2 ^-e~ kz ° (ik + z) x m is the dipole's radiation. 
The continuity of magnetic fields can be written as 

(ik + kz) x A 1 + (ik - kz) x A' r = h*(0) + h§. (9) 

where h p (r, z) = J k e lkr h p (k, z), h (r, z) = 
J k e lk r+pz h.Q. Applying (ik — kz) on M eliminates A r 
and with (Vxh) 2 (r,0) = and using ( V x h p ) z (r, 0) = 
the total field can be written as 



h(r,z) 



IkJz' A 

2 p{z+z '\p- k)(z-ik)} 



ik 



-[(ik x z) ■ A l ](kz + ipk). (10) 



Applying (ik + pz) on Q identifies, after some algebra, 
our goal 



h r (r,z) 



[S] 



A 2 



c ik-[r-u(z')] e - 



-kz+pz'j^ -fll) 

k + p 1 ; 



where [S] = J k e^ r - vz (^(ik x z) ■ A 1 ) [ik-z] stands 

for a pure superconductor response, i.e in the absence of 
a vortex; this term does not contribute to the magnetic 
fluctuations. 



The form of this Gree n's function reproduces the effects 
of image vortices^^ and allows for straightforward cal- 
culations for a general shape vortex. Corrections to ^ 
and to the following results are of order 0(e~ L / x ) for a 
finite thickness L of the superconductor. 

It is convenient to shift the particular solution of ^ 
h p (r, z) to hp = hp — A 2 V(V • hp) so that V • h p = 0, 
hence 



A 2 



du . 



az' 



+ X 2 e^ z+Z '\ik + pz)] 



(7) 



III. VORTEX EQUATION OF MOTION 

We derive here the forces on the vortex: an elastic 
force, a Lorentz force due the source field, and add fric- 
tion and pinning forces. The intrinsic forces can be de- 
rived from a London free energy in terms of Vtp (Eq. [3]) 
as a variation on u(z). The result is the familiar Lorenz 
force 



F»W = ^(Vxh)x(z + g)n 



(z) 



\z) 



(12) 
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FIG. 1: Noise components as function of zo,uq, plotted for zo/X > 0.1. Parameters are a = 0.1, k = 100. Note that $ z „ 
changes sign and $ uu has a minimum, both at Zq ~ 2uq, while <J?zz(0,tio) = <j>zu(0, Uo) = 0. 



where 1 1 denotes the x, y components and (V x h) is eval- 
uated at r = u(z). Here Ff? urce (z) is proportional to 
the source A 1 while F„" s (z) is the intrinsic elastic force 
which is independent of A 4 . F„" s (z) is found from the 1st 
term of ( 10 ) and is evaluated by linearizing in u(z') — u , 
Taylor expanding u(z) — u(z'), integrating z' and keeping 
only terms that diverge as £ — > 0, 



*(*) 



(4ttA) 2 



dz 



kdk 



d 2 u 

~dz* 



kdk 



13) 



The first term diverges at z — as l/£ which is the 
upper limit on the k integration. The use of this upper 
limit is a qualitative description and is valid only for In £ 
terms, where a change in the £ coefficient is a relatively 
small correction. To be consistent with the starting Eq. 
([2]) we must eliminate the large l/£ force and impose a 
boundary condition 



du , 

dz 



1 2=0 



0. 



(14) 



We note that that the sign of this force pushes the vor- 
tex towards (14 1, i.e. the latter is stable. This boundary 



condition has the intuitive interpretation that the current 
~ V0 — A (Eq. ^ near the vortex core is dominated 
by V0, which in turn is perpendicular to the vortex di- 
rection (Eq. [3]) . Since this current must be parallel to the 
surface at z — Eq. (14) follows. In fact this boundary 



condition has been used previously in conte xt o f He ifl^ 
and superconductors^. Our derivation of ( 14 1 is rigor- 



ous (for large k) as it includes the full A and reflections 
from the surface. For k that is not large one has to go 
beyond the phase only description of Eq. ^ and use 
coupled phase - amplitude equations^. The 2nd term 
of ( 13 ) is ~ In k and in the following we take it as the 
dominant term in the elastic force. The result is then the 
well known elastic coefficient for bending of the flux line 
due to change in its length^ 



(4ttA) : 



In k 



d 2 u 



(15) 



The external source contribution to the force is found 



from the 2nd term of ( 10 ) and in terms of the source m 
it becomes 



2 W 



a ik-uo+pz— kzo 



ik 



{{ik + z) 



(16) 



The dynamics are dominated by a friction term 
where 77 can be estimated from the Bardeen-Stephen 
resulPSl (BS), 77 = 0o/(27r£ 2 c 2 p„) where p n is the re- 
sistivity of the normal state. Finally we add a pinning 
term a[u(z) — Uq] that attempts to fix the vortex at the 
location u , i.e. a columnar defect. For this strong pin- 
ning case one can estimate a as the condensation energy 
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density, hence denning a — |a0§/(47r^A) 2 strong pinning 
has a ~ 1. The equation of motion is then 



^-ln K ^ 
(4ttA) 2 dz 2 



iurju — a(u — Uq) 



(17) 



The boundary conditions are ( 14 1 and u 
— oo. The solution is, defining F^ urce (z) = 



-> Uo at z 

f k F k ee*, 



u(z) - Uq 



I In / 



Jut? -VP ( c 



Ik 
ple z/l 



(18) 



-here £ = ^^A) 2 . 



IV. RESPONSE 



The response is identified from Eq. ( 11 ), expansion in 
u(z) — u and the solution (18). The angular integra- 
tions in Eqs. ( 11|18 1 can be done analytically, leading to 
a double integral on fc, k' that is evaluated numerically. 
The relative significance of the friction and pinning is 
controlled by the ratio 



U)T) 

a 



16ttA 2 o 
c 2 ap n 



(19) 



where the BS frictiorPSl is used. For typical parameters 
of magnetic traps and type II superconductors we use 
A R3 lOOnm, to « 1MHz and p n u) « 10~ 12 , hence we esti- 
mate the ratio (19) as 10 _5 /a. Therefore, unless pinning 



is very weak we can expand in this ratio, leading to a 
response linear in ujrj. The response has then the form 



Ima 



(i) 



A 3 ac In k 

®zz(u , Zq) 



(20) 







i(uq,z ) 





where the cartesian axes are chosen in the directions 
of z,u.q,z x Uo. We note that there are two finite off- 
diagonal terms & U z(uq, zq) = <& zu {uq,zq), that are equal 
as required by Onsager's reciprocity theorem. We show 
the 4 components of the noise ^ij(u ,z ) in fig. 1. 
We note that an m z dipole placed above a vortex (i.e. 
uq = 0) does not produce a force on the vortex hence 
®zz(0,z o ) = (vanishing as ~ u 2 ,) and $ zu (0,z ) = 
(vanishing as ~ |mo | ) ■ Furthermore, & uu (u , z ) has a 
minimum and $ zu (uo, zq) changes sign, both at Zq = 2u 2 
at large zq , uq while for small zq this continues at a some- 
what smaller zq/uq ratio. For z 2 



u 2 » A 2 we find 



$zz 

$zu = f{a) 



/(«) 
/(«) 



9(z u ) 2 



(4 

{4 



u 2 f 
2u 2 f 



i-4 

3z -» (z 2 - 2up 



4>, 



/(a) = 



/(a)- 



1 



{4 + u 2 f 

A 6 1 + 45 
2(1 + V&) 4 



(21) 



where a = an 2 / In k. Hence all <$>ij(uo, zq) decay rapidly 
beyond A, with asymptotic forms $^(tio, zq) ~ z^ a , with 
a = 8, 7, 6, 6 for the zz, zu, uu and z x u components, 
respectively. 
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FIG. 2: Dependence of noise on the pinning strength a for 
2 = Mo = 0.8A and k = 100 where a = 2a(47r£A) 2 /</>o- The 
4 curves correspond to the 4 noise components as ordered in 
the inset. 
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FIG. 3: Ratio of the two noise eigenvalues in the zo, «o plane. 
Parameters are a = 0.1, k = 100. 



For not too weak pinning, a 3> In k/k 2 , we have in Eq 
and $ 



l/< 



(21) f(a 



— > 



t^-t In k, 



1/a, hence Ima 



(i) 



2&K 2 " i ' u ' =V — ~tj 

2 . Fig. 2 shows the a dependence at an intermediate 
scale, displaying a somewhat stronger decrease with a. 

The response in the z — u plane can be diagonalized 
leading to Fig. 3 with the ratio of upper and lower eigen- 
values. We find that the magnetic noise is reduced by 
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i Minimum Response Direction 
= No force on vortex at z=0 



Z /\ 1.5 




Ug/X 



FIG. 4: Direction of the eigenvector with minimum magnetic 
fluctuations (full lines). The white insertions are the direc- 
tions for which there is no force on the vortex at z = 0. 
Parameters are a — 0.1, k, — 100. 



at least a 10 2 in the direction of the eigenvector with the 
lower eigenvalue. This direction is shown in Fig. 4 by the 
thick blue lines. We find that these directions correlate 



well with the direction of m such that the source ( 16 ) 
vanishes at z = 0. The latter directions are shown as the 
inner white lines in Fig. 4. We note that for z 2 + Uq A 2 



the lower eigenvalue of Eq. ( 21 ) actually vanishes and its 
eigenvector in the (zq, uq) plane is ~ [2u§ — z 2 , 3zoito], in 
agreement with Fig. 4. 

To appreciate the scale, we note that for a normal 
metal the response is (for z < S) Ima™? — 2zo^ P nence 
the ratio is (similar for other components) 



Ima™ „ 



128tt zo 

— j ^zz(U ,Z ) — 

a In k A 



(22) 



For z m A and a m 1 the overall magnitude, as seen in 
Fig. 1, is comparable to that of a normal metal. Hence 
at short scales zq,uq < A the vortex fluctuations lead to 
relatively low noise only in the special directions . At 
large distance zq ^> A the noise decreases rapidly as z§ a 
with cr = 6 — 8 (Eq. 21 ) for the various components, 



while the normal metal's noise decreases as 1/zrj- Hence 
in the range A <C z <C S the vortex noise is considerably 
less than that of the normal metal. 

Finally we consider the case a = corresponding to 
weak pinning a <C tjlo. This case can correspond to 
weak pinning from point defects, such that the persis- 
tence length of the vorte^P^ L c ^> L (and L S> A for 
our solution to hold). In this case the vortex elasticity 
overcomes the pinning and the vortex remains essentially 
straight in equilibrium. We define a dimensionless pa- 
rameter f^-) 2 ^ = which with the BS friction, 
k = 10 2 and typical parameters as above is f=s 10~ 2 . 
For small fjuj we find that the response has the form 



Ima 



(2) 



4>o 



$ij(u ,z ) where <&y(«o,zo) is V 



independent. In Eq. (171 77 provides a restoring force, 
hence a divergence 



1/^/77 as 77 — > 0. Therefore, for a 
given a and as 77 is reduced, in the first range a -C "qui the 
noise increases, becomes maximal at a ~ rjuj, then in the 
regime a 3> rju the noise decreases as Ima-^ ~ cur] /a 2 
and finally vanishes at 77 — > 0. 

Diagonalizing $y leads to an eigenvalue A_ = with 
an eigenvector which is very close to that of the strong 
pinning case in Fig. 4. We find that 3>y and &ij have 
a very similar Zo,Uq dependence. In particu lar for Zq + 
^ we nnd that is given by Eq. (21 1, except 
for the replacement f(a) — > A 6 . Hence, comparison with 
the strong pinning case yields 



Ima 



(2) 



amV2n 2 $ij(u ,zo) 
Ima {1) ~ (2fjw) 3 / 2 \iiK<!>i j (u ,z ) 



(77a;) 3 / 2 In k 



(23) 



where the limit corresponds to z 2 + ^> A 2 . For a not 
too small this ratio is large, in particular due to the n A 



factor, originating from Ima 



(i) 



1/a 2 



V. CONCLUSIONS 

We present here a systematic treatment for a vortex 
response to an external field. Our formalism provides 



a rigorous derivation of the boundary condition ( 14 1 for 
an arbitrary vortex shape. We apply our results to the 
problem of dephasing and spin flip in cold atom traps. 
The single vortex provides an efficient tool for trapping 
cold atomiP^HIH hence the significance of evaluating the 
fluctuation effects. 

The single vortex in our system breaks translational 
symmetry, hence the reflection from the surface is not 
specular and off-diagonal elements appear in the response 
matrix. Due to these elements, there is a special [zq, uq) 
dependent direction, for which the fluctuations are con- 
siderably reduced. 

For a trapping static field in the n direction, the mag- 
netic fluctuations ~ Ima„ measure the fluctuations in 
the energy difference of two trapped magnetic subl evels 
with i in the n direction, leading to dephasing^l. By 
choosing n parallel to that of the minimum noise, Fig. 4, 
dephasing will be considerably reduced. In contrast, the 
sf rate depends on field fluctuations perpendicular to n, 
i.e. the directions in ^ of Eq. Q. Hence by choosing 
the trap direction perpendicular to those in Fig. 4 the 
noise will be reduced by a factor sw 2. 

We find that strong pinning, e.g. as from a columnar 
defect, is significant for reducing magnetic noise. Fur- 
thermore, we find a strong decay at zq > A of the mag- 
netic fluctuations ~ with a — 6 — 8 for various 
noise components, Eq. (21). The regime A -C z <C S 



is a regime where the vortex static field B is still signifi- 
cant for trapping^, i.e. B z <~ l/z 2 ,B u ~ 1/zq. Hence in 
this regime the magnetic fluctuations are significantly re- 
duced, allowing for efficient trapping with low dephasing 
and spin flip rates. 
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